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We present a geometrical description of new canonical d-dimensional codimension one quasiperi- 
odic tilings based on generalized Fibonacci sequences. These tilings are made up of rhombi in 2d 
and rhombohedra in 3d as the usual Penrose and icosahedral tilings. Thanks to a natural indexing 
of the sites according to their local environment, we easily write down, for any approximant, the 
sites coordinates, the connectivity matrix and we compute the structure factor. 
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O , INTRODUCTION 

(N 

^~^, Quasiperiodic tilings have been widely studied over the past decades, their main interest lying in their guessed 

, ^ ' relation with real quasicrystalline alloytJ atomic structure. Indeed, most of the specific features and physical properties 
associated with quasiperiodic order (electronic structure, localized phason degree of freedom, nature and dynamics 
of dislocations, ...) are already epesent in these simplified models. The two standard ways to generate these tilings 
are the cut and project methodcia which relate them to selected pieces of higher dimensional lattices, and the grid 
algorithmQ which insists on a dual picture of intersecting lines or planes. In both cases, one obtains structures made 
fvi ■ up of different types of tile which are segments in Id, rhombi in 2d and rhombohedra in 3d. The simplest (and most 
J> famous) examples correspond to tilings with two differents tiles (the one-dimensional Fibonacci chain, the Penrose 
C^ (pentagonal)l3 and Amman-BeenkerQ (octogonal) two-dimensional tilings and the icosahedral 3d tiling) although other 

C , quasiperiodic structures with more different types of tile have also been proposed. In d > 2, all these tilings display 
vN ■ a complex local order with various coordination number. 

Concerning the study of the physical properties (spectrum, eigenstates, conductipiy,...), it has been proved more 

efficient to consider simpler systems such as direct products of periodic chaingLrll3 embedded in a quasiperiodic 

^^ ' potential. The main advantage of these structures is that their characteristics (sub-ballistic transport, self-similar 

■^ . eigenstates) can be easily deduced from the Id case ; on the other hand, since all the sites have the same coordination 

C^ ' numberj-their ttmology is trivial. In 2d, Penrose-like tilings have nevertheless been studied revealing the same type of 

H , featureslljilj'' Oill-|t>ut in 3d, the incidence of the topological quasiperiodicity on the electronic properties has been 

I ' poorly investigatedlljO essentially because of the geometrical complexity. In this context, it would be appreciable to 

''O ^ have simple structures susceptible to provide a better framework for such studies. 

^ ■ The aim of this paper is to present new quasiperiodic tilings (of arbitrary dimension) following from a natural 

j^ ' extension of the Fibonacci sequence, which allow for a coherent indexing of the sites while keeping some of the 

• • ,; interesting self-similar properties of the Penrose-like tilings. We||give the procedure to build any approximant of these 

.^ ' topologically non trivial tilings using the conurahering schemec^Eil. In the quasiperiodic limit, these tilings are closely 

K^ , related to those initially considered by Rauzycj in two dimensions and will therefore be called generalized Rauzy 

Vh ' tilings. In the first part, we recall the principle of the conumbering scheme that allows to generate all the sites of any 

5t ' codimension one structures by iterating the so-called generating vector g. As an illustration, we apply this scheme to 

build the well-known Fibonacci chain for which it is possible to exactly determine the coordinates of g. In the second 

part, we generalize these results in higher dimension and we describe, in details, the characteristics of the generalized 

Rauzy tilings in 2d and 3d. The site coordinates are explicitely obtained in terms of the generalized Fibonacci numbers. 

In the third part, we give the connectivity matrix of these canonical tilings that is relevant for tight-binding electrons 

problems. The derivation of the structure factor is given in the fourth part. Finally, we propose, in the appendix, a 

dual point of view in which the Id quasiperiodic system presented have a high codimension although based on the 

same type of sequence. 

I. THE CONUMBERING SCHEME : APPLICATION TO THE FIBONACCI CHAIN 

In the cut and project algorithm commonly used to build quasiperiodic structures, one considers a D-dimensional 
hypercubic lattice and d-dimensional subspace that define the "physical" space. The {D — d)-dimensional subspace 
define the perpendicular space whose dimension is called the codimension. Thus, d-dimensional codimension one 
tilings are generated from a (d-l- l)-dimensional hypercubic lattice. The unidimcnsional character of the perpendicular 
space of these tilings provides a natural ordering of the sites which amounts to classify them according to their local 



environment : two sites with close coordinates in the perpendicular space have a similar neighbourhood. As shown 
below, each site of a codimension one tiling can actually be indexed by a unique coordonnate, its conumber, that is 
related to its position in the perpendicular space. For any approximant structure, this is achieved by usmp-a generating 
vector which, upon simple iteration, fully determines the coordinates of the sites inside the unit cellcJO. Since this 
method is valid for any approximant, the quasiperiodic structure can therefore be approached asymptotically. 

To illustrate this mechanism, we shall focus on the 2^1 case. We consider a square lattice, and draw a line 
through the origin O, denoted E" (the parallel or physical space), of rational slope a — p/q, with p and q mutually 
prime (see figure |l|). The sites of the approximant structure are obtained by an orthogonal projection on E" of the 
square lattice sites contained inside the semi-open band generated by sliding the unit square along ij". Since i?" has 
rational slope, it contains a set of regularly spaced sites of the square lattice which defines a unit cell with n = \p\ + \q\ 
sites and a cell vector A". We define the so-called generating vector g as the vector of smailest norrrO, joining the 
origin to the square lattice site closest to E^^ (but not belonging to i?"). It can then be shownc3£2l that the sites inside 
the band are obtained by successive translations of g {modulo A") : 

1""' — j S modulo a" , j e [0, n — 1] . (1) 

This indexation by a unique coordinate (the conumber) provides a natural ordering with respect to the distance 
(before projection) from the parallel space, i. e. with respect to the local environment. Note that the origin has, by 
definition, the conumber 0. For arbitrarsi-Sjeauences of approximants, the generator coordinates are given by a specific 
element of the Farey tree decompositionE3^Eil but, for particular rational slopes, it is possible to explicitely write down 
the generating vector as we shall now see for the Fibonacci chain. 

We denote by S2=(ei, 62) the canonical orthonormal basis of the square lattice, and consider the vector A^ whose 
coordinates in B2 arc given by {Fk, Fk~i). The index k refers to the approximant order and the Fibonacci sequence 
{Fi)i£z is defined as follows : 

Fi+i =Fi+ Fi^i with Fo^Fi = l . (2) 

Usually, this sequence is only defined for / > but, here, we shall also consider negative values of /. Aa-readily seen 
on the recursion relation (0) lim Fn+i/Fn = r where the golden mean r = (f + \/5)/2 is the PisotEJ solution of 

the quadratic equation x'^ = a; + 1. Aj, is then defined as the vector of smallest normcJ with integer coordinates 

perpendicular to A^. (AJI, A^) forms a basis of the so-called trace lattice A^. One also defines the band lattice Sfc 

generated by AJI and the vector u = (1, 1) (see figure (|^) for k = 4). Let Lk (resp. 5*^) be the matrix transforming 
the canonical basis B2 into the A^ (resp. S^) basis. The number of sites of the square lattice contained in the A^ 
(resp. Sk) unit cell thus reads It — \ detL^I (resp. Sk — \ det5fc|). To determine the coordinate of the generating 
vector gfc, we introduce the matrix M : 

M^(;j), (3) 

that verifies the golden mean equation : 

M^ = M + I . (4) 




FIG. 1. Conumbering of the sites of the 4* order approximant of the Fibonacci chain upon iteration of a generating vector. 
Af = (5, 3), a| = (-3, 5) and g4 = (2, -3). 

The powers of M can be expressed in terms of the Fibonacci numbers : 

M" = ( ;," ;?-' I , Vn e Z. 



Fn-1 Fn-2 



(5) 



Note that we also ahow the index n to be negative. Since the first column of M^ represents the coordinates of A^, 
it comes immediately that the second line of M~^ gives the coordinates of AJJ, in the basis B2 '■ 



AJI =(i^_fe_i,F_ 



k-2) 



The matrices Lk and Sk thus reads 

Lk = 



F-k-i Fk 
F-k-2 Fk-i 



and Sh = 



F-k-i 1 

F-k-2 1 



(6) 



(7) 



To determine the number of sites in the S^. unit cell, one can either directly compute s^. or one can remark that 
this number is also given by the scalar product of u and A^ . One thus readily obtains : 



Sk 



— u.Aj, — Fk + Fk-i — Fk+i ■ 



(8) 



To determine gfe, one uses the fact that the quadrilateral generated by AjJ, and g^ should have a unit area since 
it cannot contain a lattice site (except on its boundary). Since AJJ, corresponds to the second line of M^^ and since 
det-M = — 1, the first line of M^^ can be chosen to be the coordinates (in the basis B2) of the generating vector : 



gfe — iF-k,F_k-i] 



(9) 



The coordinates of the sites contained in the S^ unit cell (sites of the approximant structure before projection), 
read in the basis B2 '■ 



ri = j Sk modulo a| , j G [0, Sk - 1] 



(10) 



The modulo operator allows to carry back all the sites in the same elementary cell of the band lattice. In the basis 
(AJI, A^) of the trace lattice A^, the coordinates finally write : 

r^, = Frac (j L'^ gk) , J G [0, s^ - 1] , (11) 



where Frac (r) represents the fractional part of the r coordinates. The interest of introducing the trace lattice is that 
the coordinates of the sites after projection onto E'l are readily given by the first component of r;^ expressed in the 

(aII,A^) basis. 

In the next section, we propose to generalize the above construction to codimension one structures in any dimension, 
such that the perpendicular direction is related to generalized Fibonacci numbers . These numbers are obtained from 
the Pisot solution of the polynomial equation x^ ~ J2j=o ^"'' ^^^ D = 2 case giving the golden mean. 

II. THE GENERALIZED RAUZY TILINGS 

The first geometrical cpHstruction based on the roots of this equation (for D = 3) was proposed in 1982 by the 
mathematician G. RauzyE3. Detailed analyses, oLthese original tilings focusing on their self-similar properties and 
their fractal boundaries can be found in Refs.E^ca. Nevertheless, the tilings that we describe thereafter are different 
although close to those initially studied. Indeed, the construction proposed by Rauzy does not rely on the standard 
cut and project algorithm since the sites of the cubic lattice chosen for the projection step are contained in a cylinder 
centered around the perpendicular direction and not in the band lattice. In addition, the construction of the approx- 
imant structures is completely new. 

In order to show the simplicity of the generalization, we will use, in the following, the same notations for the 
D ^ D — 1 tilings as for the 2^1 tilings. 

A. The two-dimensional case 

As for the 2^1 case, we shall use the conumbering scheme to generate the 3^2 approximant structures. We 
endow the standard cubic lattice of a canonical orthonormal basis B3 = (ei, 62, 63) and, by analogy with the Fibonacci 
chain, we choose as perpendicular space the direction defined by vector A^ = {F^jF^-i, Fk_2) where the generalized 
Fibonacci sequence {Fi)i^z is defined as follows : 

Fi+i =Fi+ Fi^i + Fi^2 with F_i = 0, Fo = Fi = 1 . (12) 

As previously, the ratio of two successive elements of this sequence converges toward an irrational limit : 

4 + (19 + 3\/33) '^' + (19 + 3733)"^' , , 

lim F„+i/F„ = a= ^ '- ^- '- ~ 1.83929 , (13) 

3(19 + 3V33)^ 
where a is the Pisot root of the cubic equation x^ ^ x^ + x + 1. We introduce the matrix M : 

(14) 




which satisfies A'P = I\'P + Af + 1 , and whose successive powers read 



Fn Fn-l+Fn-2 Fn-1 

M"=l Fn-1 F„_2+F„_3 Fn-2 I VneZ. (15) 



Fn-2 Fn^3 + Fn-4: Fn- 



3 



Below are shown the first (positive and negative) powers of M which display a remarkable pattern on both parts of 
the identity matrix M^. 
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(16) 



Note that the vertical sequences of numbers, in each column, obey relation (O) from bottom to top. To build the 



fc*'* order approximant, one remarks that A^ corresponds to the first column of M'' so that the two vectors All and 



2 



II —h 

AJl which generate the parallel space are directly obtained from the second and third line of M : 

A^ = (^_fe_i,F_fe_2+-F-fc-3,i^-fc-2) (17) 

Af = (f_fc_2,i^-fc-3 +^-fc-4,i^-fe-3)- (18) 

111 II 2 II 1 II 2 

(A)i , Af. , A^) defines a basis of the trace lattice A^ ; the band lattice S^ is generated by (A)! , Aj, ) and the vector 
u ~ (1,1,1), which joins the origin to the extremity of the unit cube whose projection onto A^ has the highest 
positive magnitude. The number of sites in a S^ unit cell is given by : 

Sk = u.A;^ = Ffc + Fk-i + Fk-2 = Fk+i. (19) 

The generating vector gfe is determined by the condition : det f AJl , AJI , gfc j =1. Since det(M^'^) = (dot M)^^ = 1, 
one can identify the coordinates of gfc with the first line of AI~'' : 

gfc = (F_,., F^k-i + F_fe_2, i^-fc-i). (20) 

and the coordinates (before projection) of the sites contained in a Efe unit cell read in the basis B3 : 

ri = J gfc modulo (Af , Af ) , j e [0, Sk - 1] . (21) 

II 1 II 2 
It is still interesting to express these coordinates in the trace lattice basis (A)i , A)i , A^) : 

ri = Frac (j L^^gk) , J e [0, s^ - 1] , (22) 

where : 

F-k-i F^k-2 Fk 

F-k-2 + F^k-3 F^k-3 + F-k-i Fk-i I . (23) 

F-k-2 F_k-3 Fk-2 

After projection onto _E[|, the j*'' site coordinates, in the basis (AJI , AJI ) are given by the two first components of 
r^,. Figure H displays an elementary cell of the 10*'' order approximant. 




FIG. 2. Unit cell of the 10*'' order approximant (504 sites). The first iterations of the generator are shown. 

In the initial cubic lattice, each site is six-fold coordinated so that, after projection, the coordination number z is 
such that 3 < z < 6. However, as discussed below, the tilings generated here only have sites with 3 < z < 5. 



B. The higher dimensional case 



The simplicity of the above construction results directly from the particular properties of the generalized 
Fibonacci numbers. It can therefore easily extended to codimension one tilings of any higher dimension. Thereafter, 
we briefly describe the three-dimensional case and the d-dimensional case. 

The construction in the 4 — + 3 case is based on the Pisot root of the equation : 



whose approximants define the perpendicular space direction in the 41? space. 
The associated generalized Fibonacci sequence (-Fi)igz is defined by : 

F,+i = F, + Fz_i + F,_2 + i^i-3 with F_i = F_2 - 0, Fo = i^i = 1 , 

and the M matrix satisfying (pj) reads : 



M 



/ 1 1 1 1 
10 
10 

Vo 1 0, 



(24) 



(25) 



(26) 



so that 



M":= 



Fn Fn-l + Fn-2 + -Fri-3 

Fn-1 Fn-2 + Fn-3 + Fn-4 

Fn-2 Fn-3 + Fn-i + -Fri-5 

-^n — 3 -^n — 4 i J^n — b i -^n — 6 



Fn^l + Fn-2 Fn-1 

Fn~2 + Fn-3 Fn-2 

Fn-3 + Fn-i Fn-3 

Fn-4 + Fn-5 fn-4 



Vne Z 



(27) 



The fc order approximant is obtained by choosing the vector A^ = (F^, Fk-i, Fk-2, -Ffe-s) whose coordinates are 
given by the first column of A/*^ . The three parallel space vectors reads : 



A^ = (F_fc_i,F_fc_2 +F-fc-3 + -F-fc-4,F-fc-2 +F_fe_3,F-fe-2) (28) 

Af = (F_fe_2, F-fc-3 + -F-fe-4 + F_fe_5, i^-fe-3 + F_fe_4, F_fe_3) (29) 

Af = (^-fc-3, J^-fc-4 + F-fc-5 + F-fc-6, i^-fc-4 + f-fc-5, i^-fe-4) , (30) 

111 II 2 II 3 

that are respectively the second, third and fourth hne of M . The band lattice Efc is generated by {A'j, , Ajl , A)l ) 
and the vector u=(l,l,l,l). The number of sites in the E^ unit cell (and therefore in the approxiniant unit cell) is 
given by : 

Sk = u.A^ = Fk + Fk-i + Fk-2 + Fk-3 = Fk+i, (31) 

and the coordinates of the generating vector are given by the first line of M^'' : 

gfc = {F-k,F-k-i + F-k~2 + F-k-3,F-k-i + F-k-2,F^k-i) , (32) 



The sites coordinates have an expression similar to Eq.(p2[) where L^. is, as before, built from the coordinates of 

Ai\Af,Af andA^. 

The generalization in any dimension D is straightforward considering the Pisot root of the equation : 

D-l 

x^ = ^ x^ , (33) 

The remaining of the construction follows mechanically. 

Let us end this section by remarking that nothing prevents, in any dimension, to map the whole structure along a 
direction such that all tiles (rhombus in 2d, rhombohedra in 3d) are identical up to rotations. Indeed, once the sites 
have been generated with a given choice of the perpendicular space, this can be achieved by projecting perpendicularly 
to the direction of the vector u whose each coordinate equals 1. In that case, all the sites are mapped onto the full 
set of simple lattice vertices (triangular in 2d and centered cubic in 3d) with missing edges. We are then left with 
a quasiperiodic decoration of a periodic structure, which has the same connectivity but a different geometry as the 
above discussed tilings. Such deformations might be useful in some context, like for instance when one aims to build 
row-by-row or plane-by-plane transfer matrices, in which case the preexisting reticular stratification of the lattice can 
be very helpful. 

III. CONNECTIVITY MATRIX OF CODIMENSION ONE TILINGS 

In this section we write down the so-called connectivity or adjacency matrix K of codimension one tilings. This 
matrix is defined by the following rule : Kij = 1 if sites i and j are nearest neighbours {i. e. connected by an edge) 
and zero otherwise. Such a matrix is of special interest for the study of the electronic or phononic excitations in 
tight-binding like approaches. In this case, the hamiltonian is, indeed, simply proportional to K. Being able to easily 
obtain such a matrix may therefore be highly valuable. 

The idea consists in labelling the sites with respect to their conumber. As discussed in section II, this indexation 
classifies the sites according to their coordinate in the perpendicular space direction, or equivalently, in terms of 
their local environment. This implies that the conumber difference between two nearest neighbour sites along a 
given direction is a constant modulo the number of sites in the unit cell. This difference is readily obtained from 
the projection onto the perpendicular space of each edges meeting at any site of the initial hypercubic lattice. It 
is therefore sufhcient to determine, for a given site, the conumber of its nearest neighbours, to write down the full 
connectivity matrix as a band (Toeplitz-like) matrix. 

We shall first illustrate this method with a Fibonacci chain approxiniant, defined by the perpendicular space vector 
A^ = (Ffe, Fk-i) and determine the conumber of the nearest neighbours of the origin. In the canonical basis B2 of the 
square lattice, these sites have the following coordinates (0, ±1) and (±1, 0). Their projection onto the perpendicular 
space is given (up to a sign), in the trace lattice A^ basis, by A-^/lk- In addition, and by construction, one has 
A^.gfe = 1/lk in the A^ basis. One therefore deduces that the conumber of the sites which are nearest neighbours of 
the origin (once carried back in a unique unit cell) reads Fk et Fk-i- 



As an example, we give the connectivity matrix K of the 4*^ order approximant of the Fibonacci chain displayed 
in figure (|l|) with periodic boundary conditions: 

/O 1 1 0\ 



K = 



10 10 


10 1 


10 10 


10 1 


10 10 


10 10 


10 10 0/ 



(34) 



There are indeed F^ = 8 sites in the elementary cell and the conumber of the origin neighbours are F/^ — 5 and F^ — 3 
(see fig. p. 

The above discussion applies to the generalized Rauzy tilings : their connectivity matrix remains a band matrix 
once the sites have been ordered according to their conumber. The conumbers of the origin's nearest neighbours are 
still given by the coordinates of A^ which are expressed in terms of the generalized Fibonacci numbers. We display 



below the connectivity matrix of the 4*^* 
periodic boundary conditions : 



order approximant of the two-dimensional generalized Rauzy tiling, with 



K 



/O 1 1 
10 10 
10 10 1 
10 10 
10 10 1 
10 10 
10 10 
10 10 10 
10 10 1 
10 10 
10 1 
10 

\0 1 



1 0\ 

10 

10 

10 10 
10 1 







10 10 1 
10 10 
10 10 
10 1 
10 10 
10 1 
10 10 
10 10 0/ 



(35) 



One has F5 = 13 sites in the unit cell, and the conumber of the origin's nearest neighbours are F4 — 7, F3 = 4, and 
F2 = 2. One can therefore appreciate the advantage of this family of tilings that provides topologically non trivial 
structures which are nevertheless very easy to construct and to encode. Especially in 3d, this is undoubtlessly one of 
the simplest type of connectivity matrix that one could expect for quasiperiodic structure approximant of any order. 



The last point we would like to discuss concerns the proportion of the different types of site and their coordination 
number. Clearly, the coordination number of site j is the number of 1 in the j*'* line of the connectivity matrix K. 
Hence, in the above example one finds four 3-fold coordinated sites , five 4-fold coordinated sites, and four 5-fold 
coordinated sites but no 6-fold coordinated sites. It is easy to show that the condition to have sites of maximal 
coordination number z = 2D in a generalized Rauzy tiling of type D ~^ {D ~ 1) reads : 



Ffc+i - 2Fk > ^ 



Fk, 



>2 . 



(36) 



But looking to the definition of the generalized Fibonacci sequence {Fi)i^z in the _D-dimensional case 



D-l 



n+l 



Y^ Fn-j with Fj = for j G [2 - D, -1] and Fo = i^i = 1 



(37) 



one can easily prove, by recursion, that the condition ( P6|) is rpyer fulfilled. This is in contrast with the two-dimensional 
tilings originally introduced by Rauzy and discussed in RcfB3 that displayed 6-fold coordinated sites. 



IV. STRUCTURE FACTOR 



The simple form taken by the codimension 1 tiUng sites coordinates is very helpful to compute the structure factor 
(which amounts here to the Fourier transform of the siter-distribution) . Indeed, the existence of a generating vector 
allows to write the structure factor as a geometric seriecl. Denoting by n the number of sites in a unit cell and 
assuming that the sites (atoms) have the same form factor, the structure factor S reads : 



_. n— 1 



(38) 



where q is a reciprocal space vector of the approximant structureR and where r-' represents the vector associated to 
the site j (after projection). In the codimension one case, we have seen that the coordinates were simply expressed 
in terms of the generator g : 



r^ = jg, i e [0,n- 1] 



(39) 



This expression slighly differs from (|2l| ) since here, r^ and g denotes the vectors after projection. In addition, the 
sites are not carried back in the unit cell through the modulo operation, but since one only considers reciprocal space 
vector, this does not affect the result. Thus, the structure factor S simply writes : 



n-l 






i=o 



The peak intensities (think of a diffraction experiment) are proportional to : 



|5(q)| = - 
n 



sin (nq.g/2) 



sin (q.g/2) 



(40) 



(41) 



In the quasiperiodic limit, the reciprocal space unit cell shrinks toward zero, and the Fourier spectrum becomes 
dense. However, going from an approximant to the next one mostly amounts to add new peaks of smaller and smaller 
intensities while letting almost unchanged the previous peaks. This is why any threshold function filtering out the 
peaks below a certain (arbitrary) value gives rise to a point-like diffraction pattern (see figure (0) for an approximant 
of a two-dimensional generalized Rauzy tiling) . 




^ Since approximant structures are periodic, only those vectors contribute to the structure factor. 



FIG. 3. Structure factor for the 10*'' order approximant (504 sites per unit cell). Only intensities greater than an arbitrary 
threshold are shown. The radius of the circles is proportional to 151. 

Note that the Fourier spectrum does not display a particular point-like rotational symmetry (except the trivial 
parity q <-> — q), in constrast with the more commonly studied octagonal and Penrose quasiperiodic tilings which 
have higher codimension (2 and 3 respectively). 

V. CONCLUSION 



The generalized Rauzy tilings form a canonical set of codimension one quasiperiodic structures. Their construction 
is based on the generalized Fibonacci sequence, whose properties allow for a rather simple method to get the sites 
coordinates as well as their coordination number. We have shown how to get closed formulas for these quantities valid 
for any approximant structures, and which therefore can be carried up to the quasiperiodic limit. The unidimensional 
nature of the perpendicular space allows one to write down the connectivity matrix in terms of Tccplitz band-like 
matrix, the position of the non vanishing band being directly determined from the generalized Fibonacci numbers. . 

These properties should rank this family of tilings among the most interesting to be studied. As shown recentlyaEEl, 
they are, to many respect, much simpler than the celebrated Penrose-like tilings, while displaying the same kind of 
physical properties. 

We would like to thank P. Arnoux, N. Destainville and L. Vuillon for very fruitful discussions about Rauzy original 
works. 



APPENDIX A: ONE-DIMENSIONAL GENERALIZED FIBONACCI CHAINS OF HIGH CODIMENSION 



Up to now, we have described d-dimensional codimension one structures with a perpendicular space built from the 
generalized Fibonacci series. But nothing prevents to build, along the same line, one-dimensional structures made 
up of {d — 1) different edges, by switching the respecting role of the perpendicular and parallel space. Let us briefly 
describe the obtained sequences and start by discussing the 3 — > 1 case. 

We write down the 3-letter substitution which is given by the matrix M, 

A^ AB, B ^ AC, C ^ A, 
which generates (starting from A as an example), the following set of sequences Si: 

Si ^A,S2= AB, S3 = ABAC, S4 = ABAC ABA, S5 = ABACABAABACAB, ... 

Note that the length of the Z*'* sequence is given by the generalized Fibonacci numbers Fj. The sequences follow 
a concatenation rule 5";+! = SiSi-iSi^2 and the number of occurences of A (resp. B, C) in the sequence Si is Fi-i 
(resp. F;_2, Fis). Note that a related sequence has actually been considered in a previous workEd but with the 
reverse concatenation rule : Si+i — Si^2Si^iSi. 

The generalization for higher codimension is straightforward. 
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